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Conductivity of the defectless, perfect crystal graphene is found at the neutrality point at zero 
temperature and in the limit of large dielectric constant of the substrate. The steady state of the 
graphene with weak current is assumed to be an ideal, rare plasma of particle and hole excitations 
governed by the Boltzmann kinetic equation. 
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The conductivity of a single-atom graphene layer as 
a function of carrier doping concentration shows a pro- 
nounced minimum at the neutrality, compensation point 
[3, 0, Here both particle and hole excitations are 
present in the vicinity of two special points in the Bril- 
louin zone where the dispersion becomes almost rela- 
tivistic, massless, cone-shaped one 0, [Bj]. In this case 
the theory predicts a universal conductivity for the non- 
interacting particles moving in a static disorder potential 
Q. However, in graphene the experimental conductivity 
exceeds by few times the theoretical universal conductiv- 
ity. Recent efforts Q to clear this discrepancy mainly 
focus on the effects of disorder and have successfully ex- 
plained a linear dependence of the conductivity as func- 
tion of the doping concentration Q]. Yet one point is 
being missed. Namely, the current of mobile charge ex- 
citations in a defectless, perfect crystal graphene can be 
changed, relaxed in the course of their mutual Coulomb 
interaction. And a question what is the conductivity of 
the defectless graphene is meaningful. Unlike the usual 
electron liquid where the current is synonymous in many 
situations to the momentum and, thus, where there is 
a rigorous current conservation by the Coulomb inter- 
action. In this paper, the conductivity of the defectless 
graphene is found at the compensation point and in the 
limit of weak Coulomb interaction constant. 

In the tight binding model on the honeycomb lattice, 
that may represent the band structure of the graphene, 
the current operator of few electrons: i XX^/xV'x+a — 
Vx+a^x), does not commute with the Coulomb interac- 
tion operator. However, close to the cone apexes in the 
momentum space, or equivalently in the long wavelength 
limit, the current conservation is approximately restored. 
Nevertheless, this symmetry can not be projected onto 
the graphene directly as the latter has two types of excita- 
tions: particles and holes, in two relevant crystal bands. 
One band is completely filled and the other is empty, 
at the compensation point, and these two bands touch 
each other at the two cone apexes. Thus, the current 
operator has two differently ordered particle and hole 



terms: j =y^^(p)de + /<9p^ + (p)-^_(p)de_/<9pV>l(p), 
as well as the inter-band current known as the Zitterbe- 
wegung. Here, p counts momenta in the Brillouin zone. 
To show how a commutation relationship changes when 
one half of the states is being filled, consider two oper- 
ators: A — Tp + Aip — tjjf Aijvbj and B = ip + Bip, acting 
in the Fock space of electron states numerated by i. If 
the matrices A and B do commute: [A, £?]_ — then the 
operators also commute: [A, B]_ = ip + [A, B]_ip = 0. 
However, if a subset of the states {i} is filled then elec- 
tron operators have to be normally ordered to represent 
excitations. In general, the ordered operators no longer 
commute: [: A :, : B :]_ = Tr ((AN B - BNA)(\ ~ N)), 
where N is the diagonal occupation matrix with the en- 
tries one for the filled subset {i} and zero otherwise. 

The main precondition for the non-conservation of the 
current is being close to the compensation point. A neu- 
tral cloud of non-interacting particles and holes responds 
to the electric field by separation: particles to one side 
whereas holes to the opposite side. In the momentum 
space, though, they all move in one direction. On the 
other hand, if a neutral particle-hole cloud is coupled 
by strong Coulomb forces it behave like a collection of 
pairwise neutral atoms. The response of these to the 
electric field is, initially, a polarization rather than a cur- 
rent. Therefore, a precise value of the conductivity is 
determined by the mutual Coulomb interactions. A mi- 
croscopic process that changes the current is shown in 
the Fig.l. A particle 1 and a hole 2, considered as a pair, 
has zero net momentum and non-zero net current. The 
velocity of an electron in the state 2 is opposite to the 
velocity of an electron in the state 1. But a hole is the 
absence of an electron, therefore, the total current of the 
pair (1,2) is non zero. In the process of Coulomb interac- 
tion the pair (1,2) can be scattered into the new position 
(3,4) with the same total momentum and energy. We 
observe that the net current of this pair in the new state 
(3,4) has been reversed. Last argument is that although 
the kinetic energy in the graphene is artificially Lorentz 
invariant the Coulomb interaction is instantaneous and 
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FIG. 1: A microscopic scattering process of the particle-hole 
pair (1,2) into the particle-hole pair (3,4) that conserves the 
momentum and the energy but changes the current. The ex- 
citation velocities are shown by the arrows. X-axis represent 
the momentum, whereas y-axis is the quasiparticle energy. 



does violate this Lorentz invariance. 

We study infinitely large, perfect 2D graphene layer 
on top of a dielectric substrate at zero temperature. Ap- 
plication of an electric field will inevitably create parti- 
cle and hole excitation due to the Schwinger mechanism. 
The work of the electric field on these excitations will 
produce the Joule heat that will be transformed into the 
lattice vibrations near the graphene layer and will even- 
tually escape into the bulk. In the balance, a steady 
distribution of particle and hole excitations will be es- 
tablished. We assume this state of the graphene to be an 
ideal, rare plasma with the excitation distribution given 
by the Fermi-Dirac function for some effective tempera- 
ture T* . The Hamiltonian in the long wavelength limit 
consists of the crystal band part [Bj]: 



a p x + a y p y , 



(1) 



where a x , a y , t z are the Pauli matrices [the first two act 
in the representation space of the crystal point group 
whereas the last one acts in the valley space] and the 
Coulomb interaction part, c = 1.1 * 10 8 cm/s 5] is the 
characteristic band velocity of the graphene that deter- 
mines the cone angle. The one-particle Hamiltonian ([1]) 
can be diagonalize by the unitary transformation: 



U = 



1 



V2 



exp it a 



(2) 



into the two crystal bands [two halves of the cone]: 
eTcr(p) = ra\p\, where p = (p x ,Py) and r, a = ±1 
are eigenvalues. At the compensation point the elec- 
tronic state of graphene is determined by a dimensionless 
Coulomb coupling: 



e 

kTic 



(3) 



where n is the half of the dielectric constants of the sub- 
strate and the vacuum. For the graphene on top of Si 
substrate g w 0.35, whereas for the graphene on top of 
Si02 substrate g 0.8. Coulomb interaction modify the 
crystal dispersion of the quasiparticle excitations [5j: 



<P) = c|p| 1 



5, Q 
4 l0g W 



(4) 



in the long wavelength limit |p| <C Q, where Q is the Bril- 
louin zone size, and this non-linearity of the dispersion 
will be important below. 

In addition to the true non-equilibrium state of the 
graphene with the current described by the electron dis- 
tribution function in the momentum space: F a (p), where 
a = ±1 specifies the two crystal bands, upper/lower 
halves of the cone, we consider also an imaginary 'equilib- 
rium' state with relaxed, zero current but with the same 
excitation energy. We do not consider here the graphene 
states that have a particle-hole coherence of any kind [see 
e.g. Q], as this may lead to the time dependence of the 
coherence order parameter averaging out its effect. In our 
solution F a (p) does not depend on either spin or valley 
indices. The spin and valley spaces give the total degen- 
eracy of electron states in the graphene N = 4. As the 
total momentum of the scattered electrons is conserved 
in crystal [neglecting the Umklapp processes], we search 
for the graphene state with zero total momentum. In this 
state there are on average as many holes as particles in 
every small momentum cell. Therefore, we can describe 
the particles by the distribution function in the momen- 
tum space F+(jp) whereas the distribution function for 
the holes 1 — F_(p) has to be the same. Thus, the dis- 
tribution function possesses the particle-hole symmetry: 
1 — F a (p) = F- a (p). In the 'equilibrium' state of the 
graphene the electron distribution function: 
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exp(a|p|/(p» + 1 



(5) 



makes the collision integral due to the Coulomb interac- 
tion to vanish for any scale parameter in the momentum 
space: (p). This scale also defines the effective temper- 
ature of the electrons in the 'equilibrium' state of the 
graphene: T* = c(p). The distribution function Eq.([5]) 
satisfies the electron- hole symmetry: f- a (p) — 1— /<*(p)- 
From now on we will isotropically rescale the momentum 
space in the vicinity of the cone points to set (p) = 1. 

The Boltzmann kinetic equation defines the steady dis- 
tribution function in the state p balancing the two pro- 
cesses - the drift of excitations in the electric field and 
the redistribution of excitations during their collision: 



eE^ = St a (p) 
op 



(6) 



In the lowest order of the Coulomb coupling g, the sec- 
ond order Fermi golden rule, the collision integral reads 
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[for short notations a 3 , a 4 of the out-going electrons are 
inverted] : 



St ai (pi) 



a2CK3C*4 



Tr r (|lCa 2 4 (PiP2P3P4)| ; 



(2tt)S(J2 o*c(Pi)) II (Pi) - II *"-«« (Pi 



(27r) 2 !5(pi + p 2 - p 3 - P4) 



rf 2 p 2 cPp 3 ri 2 p 4 
(2tt) 6 



(7) 



Below we use interchangeably the notation: pi = p, p 2 = 
p', P3 = p' + q and P4 = p' — q. The Coulomb matrix 
element is weakly screened, as the plasma of excitations 
is assumed to be rare in the limit of small (p) [with the 
Debye screening radius being large Rd ~ h c 2 /e 2 T*],: 



^ Q3 ^ 4 (P1P2P3P4) 



2^e 2 



1 



1 

2 V K |Pi - P3I 2 

2?re 2 1 - z 2 z% 1 



Z\Z% 1 — Z 2 Z. 



2 

z\z\ 



K|P2 - P3| 



<5 Tl 7- 3 8 T 



^TiT4^r 2 T3 



(8) 



where the notation = a^p* + ipi)/\pi\ is used, r in- 
dices run over the total N — 4 spin-valley degeneracy 
space of the graphene. The square of the matrix clement 
Eq.® includes two terms - the direct and exchange ones. 
The exchange term vanishes when two scattering excita- 
tions have different spins or valleys. 

In the graphene state with current the electron distri- 
bution function can be written as: 



exp 



a\p\+a(eE-p) x(|p|)/|p| 



1 



(9) 



where x(p) defines the perturbation of the distribution 
in electric field. It is better to satisfy the condition 
X(p) — > as |p| — > 0. Also, x(p) in Eq.® explicitely 
conserves the number of electrons, their total energy and 
their total momentum. We linearize the Boltzmann ki- 
netic equation Eq.Q with respect to x(p). The lin- 
earized collision integral becomes a matrix, which is sym- 
metric due to a detailed balance of the direct and time- 
reversed processes in the steady state. The current in the 
linear responce reads [the band velocity v a — ap/\p\]: 



j[x\ 



4^p(p-E)up)f- a (p) d2p 



(27T) 2 

(10) 

For strictly linear dispersion and for the collinear ori- 
entation of all momenta p||p'||q the argument of the 
energy delta- function in Eq.([7]) becomes degenerate, i.e 
ai|p| + a 2 |p'| + a 3 |p + q| + a 4 |p' - q| = for any 
|p'| and |q| provided three conditions are met: a\ = 
a 2 sgn(p-p 2 ) = -a 3 sgn(p-p 3 ) = -a 4 sgn(p • p 4 ). Taking 



the direction of vector p as x and expanding around the 
collinear configuration of four momenta we find: 



/2 



1y (Py - Ivf 
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2p> 2(p + q) 2(p'-q) 4^ 



Pi log \pi\ (11) 



where pi — pi X and where the effect of the renormaliza- 
tion of the velocity 5c = (<?/4) log(Q/(p)) is omitted. The 
last term in Eq. lfTTj) can be approximated as ±g. The 
integration of the energy delta- function 5(AE) with re- 
spect to the y momentum components gives the Jacobian 
y / pp'(p + q)(p' - q)/\p\, provided pp'(p + q)(p' - q) > 0, 
and the large logarithm 21og(l/g). Thus, in this leading 
large logarithm approximation: 



h log G 



> 1 



(12) 



the linearized Boltzmann kinetic equations reads: 



-Hoc 



X(\P\) + X(\P'\) - x(\p + g|) - X(\P' - l\) 
00 (eP + l)(eP' + l)(e-P ~ 1 + l){e~P' + 1 + 1) 

-\P\ 

{eP + \){e~P + I) 



\JPP'{P + q)(p' - q) dp'dq 



2n 



{13) 



where A = 2Ng 2 log(l/g) is the Coulomb integral, and 
the condition pp'(p + q){p' — q) > is enforced in the in- 
tegrand. The exchange term vanishes in the leading log- 
arithm approximation. The Debye screening mass makes 
the integral in Eq. (|13p to converge as the principle value 
in the vicinity of q — 0. Due to few symmetries of the 
integral in Eq . (fTB"]) : (p <-> p' , q «-> —q), p <-> — p — q 
and p' <-> — p' + q, the Eq. (|13p is a symmetric operator. 
Thus, the Boltzmann kinetic equation is the variation of 
the functional: 1Z[x] — ^[xL where 



+ OO 



\Jpp"\P + q)(p' - q) dpdp'dq 
q 2 (2tt) 2 

(X(\p\)+x(\p'\)-x(\p + q\)-x(\p'-q\)) 2 
(eP + l)(eP' + l)(e-P ~1 + l)(e~P' + 1 + 1) 

PX{p) dp 



(eP + l)(e~P + 1)2tt 



(14) 



The existence of this positively defined functional 7t[x] 
proves that the conductivity is positive. Indeed, in the 
minimum: K[x] - £[x] < because K[0] - E[0] = 0. The 
conductivity a — S[x] > 72. [x] > 0. 

The equation (|13p is contradictory and has no solution 
as the integral over all p applied to the left hand side 
is zero. It means that the leading, large logarithm ap- 
proximation is insufficient. However, from the all next 
order terms of the Boltzmann kinetic equation we need 
only their combined action on the homogeneous, in the 
momentum space, mode x(p) — Xo — const, which is be- 
ing neglected by the leading term Eq. lfTB")) . The 'leading' 
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order of this 'subleading' term is g 2 without the large 
logarithm. The mode xo arises in the process of parallel 
shift of all momenta |p + a| = p + (p • a)/|p|. We write 
the 'subleading' term as a projection: iV<7 2 |<I>(p))(<I>(p)|, 
onto a function $(p). The function <i>(p) can be found 
in the closed form as an integral. We parameterize the 
four momenta of two scattering electrons by pi — a, |pi|. 
These also define the mutual angles of the momenta upto 
to the four-fold discreet flip transformations. The colli- 
sion integral gives the function <I>(p) = St a (p)[x = 1], 
independent of a, and we set a = +1. Its exchange part 
is rather complicated whereas the direct part reads: 



Mp) = // ■ 

J J — oo 



(eV + l)(eP2 + l)(eP3 + l)( e P4 + 1) 

U 



us 



(2Q 2 



u — s) Arcth 



where p4 — —p ~p2~Pz due to the energy conservation, 
and the parameterization: Q = p + pz = —p2 — P4, u = 
ippz, s = 4p 2 p4 [with Q 2 > |u|,|s|] is used, satisfying 
the condition us > in the Eq. (fT5]) . As long as Q(p) 
is known and the homogeneous mode xo is singled out: 
X(p) = Xo+Xi(p)j the Boltzmann kinetic equation reads: 



m? 2 $(p)xo + -Hxi] 

OX 



p 



(eP + l)(e-P + l) 



(16) 



In the large logarithm limit Eq.([T2]), the function Xi(p) 
is relatively small |xi(p)l ^ lxo|- Integrating Eq. (fT6)) 
with respect to p eliminates all non-homogeneous modes 
in the kinetic equation that determines xi an d leaves the 
equation for xo only: 



Ng 2 C X o = 



-\p\ dp 



o (eP + l)(e-P + l) 



= -log(2) (17) 



Solution is: xo = ~ l°g(2)/iVCg 2 , where C is the average 
of the collision integral over the homogeneous mode: 



C = 



<&(p)dp 



(18) 



We estimate numerically the direct Cd ~ 0.69 and the ex- 
change C ex » —0.1/N parts of C = Cd + C ex . Thus, the 
distribution of excitations with current is approximately 
the same as without current just translated in parallel in 
the momentum space by a vector proportional to the elec- 
tric field. This solution, though, does not vanish in the 
momentum origin. We estimate numerically that adding 
some gap to the collision integral in the origin AS(p)S(p') 
gives a necessary crossover of x(p) to zero on the scale 
|p| < g, and as g — > this feature has negligible effect. 



The conductivity is determined by the current Eq. (jTUJ) : 
j[x] = ~ AHog(2)(e 2 /?i) xo E, neglecting small xi con- 
tribution. From this equation, we finally obtain the con- 
ductivity of the defectless, perfect crystal graphcnc: 



e 2 log 2 (2)/C 
h 2irg 2 



(19) 



in the limits of large dielectric constant of the substrate 
g^O and, also, the large logarithm log(l/<?)/(27r) — > oo. 
Remarkably, the conductivity Eq. ([T9| does not depend 
on the 'heating' scale (p) = T*/c. Inclusion of spe- 
cific mechanisms of energy relaxation due to say electron- 
phonon interaction is necessary to determine this scale. 
At any rate, T* is small in the limit of weak electric 
field T* ~ E 1 . However, the power 7 can be rather 
small in the realistic electron-phonon models. For the 
Coulomb coupling g « 0.35 the minimum conductivity 
in the Eq.(|19p corresponds to the experimental values 
around p max ~ 4 kOhm. Our numerical estimation of the 
different parts of the Boltzmann kinetic equation shows 
that the large logarithm approximation begins to work at 
around g < 0.2 whereas for the experimental condition 
g ~ 0.4 an increase of the conductivity Eq. (fT9]) by 30% 
or so has to be expected. 

As the gate voltage breaks the particle-hole symmetry 
and creates a net charge on the graphene: e(Nh — N e ), 
the electric field induces non-conservation of the total 
momentum: dP/dt = e(Nh — N e )E. This runaway evo- 
lution of the excitation distribution can not be controlled 
by their mutual Coulomb interactions alone, as they con- 
serve the momentum, and some defects violating the 
translational symmetry is required to stabilize the steady 
state. This situation can be treated in a similar fashion. 
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